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What  are  the  theoretical  and  experimental  physical  differences  between 
waves  and  turbulence?  The  motivation  behind  this  question  is  related  to  the 
practical  problems  associated  with  laser  beam  propagation  and  pollution 
transport  in  the  atmosphere.  Because  turbulence  causes  mixing  and  waves 
do  not,  one  must  not  regard  turbulence  as  a  field  of  random  waves.  The 
power  density  spectrum  of  velocity  fluctuations,  when  taken  alone,  cannot  be 
used  to  distinguish  between  waves  and  turbulence.  Its  physical  interpretation 
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2  0.  Abstract  (Contd) 

can  differ  radically  depending  upon  which  type  of  motion  is  involved.  Two 
approaches  are  used  here  to  differentiate  theoretically  and  experimentally 
between  waves  and  turbulence.  The  first  involves  the  degree  of  interaction 
between  modes.  Theoretically,  this  leads  to  a  new  interpretation  of  the 
buoyancy  length,  kg  ,  and  it  raises  questions  of  the  nature  of  the  buoyancy 
subrange  of  turbulence.  Experimentally,  it  leads  to  new  suggestions  for 
distinguishing  between  wave  and  turbulent  motion,  for  example,  by  means  of 
bi-spectral  coherence.  The  second  approach  depends  on  the  mixing  property 
of  turbulence v— Theoretically,  the  mixing  is  related  to  strong  mode -interaction, 
which  is  physically  due  to  a  vortex-stretching  cascade  process.  Experimen¬ 
tally  one  can  use  mixing  to  distinguish  between  waves  and  turbulence  by  means 
of  a  procedure  suggested  by  Busch.  The  latter  procedure  uses  the  cross¬ 
spectrum  between  vertical  velocity  and  potential  temperature  fluctuations.  It 
is  described  here  in  some  detail  and  is  generalized  for  possible  application  to 
the  case  of  two-dimensional  turbulence.  In  addition,  a  wave  cascade  model 
is  proposed  in  order  to  explain  power  spectra  observed  in  project  HICAT, 
and  a  test  is  suggested  for  this  theory. 
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On  the  Difference  Between  Waves 
and  Turbulence  in  a  Stratified  Fluid 


I  INTRODUCTION 

Velocity  fluctuations  in  the  stratosphere  that  were  presumed  to  be  turbulence 
were  measured  (Crooks  et  al1)  by  means  of  an  instrumented  U-2  aircraft  in 
project  HIC'AT.  Figure  1  shows  some  power  spectra  of  such  fluctuations  plotted 
on  a  log-log  scale.  Typical  of  all  such  uata  in  their  report,  these  spectral  display 
a  -  A  /  3  slope.  This  happens  to  be  the  same  slope  as  that  predicted  by  Kolmogorov2 
for  the  case  of  "inertial  range"  turbulence.  Some  authors  such  as  Lilly  et  al°  have 
indeed  interpreted  this  spectrum  as  being  in  the  inertial  range  of  turbulence. 

Other  authors  have  noted  that  the  buoyancy  length  is  too  small  to  permit  an 
..icrtiai  range  interpretation  (Zimmerman  and  Loving  )  and  that  a  more  likely 
explanation  would  be  that  it  is  some  sort  of  buoyancy  subrange  turbulence.  In 
other  words,  the  inertial  range  assumption  appears  to  be  contradicted  by  the  fact 
that  the  5/3  spectrum  is  almost  entirely  at  scales  larger  than  the  buoyancy  length 
(in  Figure  l).  Following  this  suggestion,  both  Weinstock'1  and  Dewan(’,  using 
different  arguments,  attempted  to  show  that  it  is  theoretically  possible  to  have  a 
buoyancy  subrange  with  a  -5/3  dependence.  This  is  in  contrast  to  the  original 
theorie  s  of  the  buoyancy  subrange  proposed  by  Bolgiano,  ‘  Shur,^  and  Lumley. 

(Keceived  for  publication  2d  October  li'82) 

Because  of  the  large  number  of  references  cited  above,  they  will  not  be  listed  here. 
See  Hcferences,  page  27. 
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INVERSE  WAVELENGTH  (cycles/ft) 

Figure  1.  Power  Spectra  From  Project  H1CAT 
(Crooks  et  al*) 


These  original  theories  and  others  are  summarized  best  by  Lin  et  al,  1  who  also 
generalized  the  theory  of  buoyancy  subrange  turbulence.  The  original  theories 
predicted  a  much  steeper  slope  for  this  subrange,  that  is,  -3  for  the  Lumley-Shur 
theory  and  -ll/5  for  that  of  Bolgiano. 

While  an  inertial  range  approach  toward  the  understanding  of  Figure  1 
appears  reasonable  at  first  sight,  it  seems  untenable  if  one  accepts  the  experi¬ 
mental  evidence  (see  Barat  and  Aimedieu,  *  Crane,  ^  Anderson,  ^  and 


Lin,  J.  T.  ,  Panchev,  S.  ,  and  Cermak,  J.  ( 1 0 fi .“ )  A  modified  hypothesis  on 
turbulence  spectra  in  the  buoyancy  subrange  of  stably  stratified  shear  flow. 
Radio  Sci.  4:1333-1337. 

Barat,  J.  ,  and  Aimedieu,  P.  (1081)  The  external  scale  of  clear  air  turbulence 
derived  from  the  vertical  ozone  profile:  application  to  vertical  transport 
measurement,  J.  Appl,  Met.  20:275-280. 

Crane,  R.K.  (1080)  Radar  measurements  of  wind  at  Kwajalcin,  Radio  Sci. 
15:383-304. 

Anderson,  A.D.  (1 057 )  Free-air  turbulence,  J.  Meteorol.  14:477-404. 
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Kosenberg  and  Dewan^')  concerning  the  average  thickness  of  turbulent  layers  in 
the  stratosphere.  These  cited  references  agree  that  the  thickness,  which  presum¬ 
ably  represents  an  "outer  length,  "  is  of  the  order  of  100  m,  which  is  several 
orders  of  magnitude  smaller  than  the  largest  scales  of  Figure  1.  While  one  piece 
of  experimental  evidence  seems  to  indicate  that  1.  u-km  thick  layers  are  not 
uncommon  (Crooks  et  al1"),  it  has  been  argued  elsewhere  (Dewan^'l  that  the 
latter  evidence  was  based  on  a  dubious  assumption,  namely  that  turbulence  was 
responsible  for  the  fluctuations  in  question.  Instead,  these  fluctuations  seem 
more  likely  to  be  due  to  waves.  This  is  supported  by  both  outer-length  consider¬ 
ations  and  buoyancy-length  considerations.  In  any  case,  one  purpose  of  the  present 
report  is  to  explain  figure  1,  with  the  aid  of  the  hypothesis  that  the  data  of 
Crooks  et  a'.  '  ’  11  represent  primarily  waves  and  not  im  rtial  range  turbulence. 
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2.  “RANDOM  WAVES”  VS  "TURBULENCE” 


Phillips considered  the  case  of  a  field  of  random  internal  gravity  waves. 
In  particular,  he  pointed  out  that  in  this  situation  the  transition  to  turbulence 
would  not  follow  the  usual  sequence  of  "instability"  (for  example,  Kelvin- 
Helmholtz  or  Rayleigh-Taylor)  followed  by  a  sudden  "turbulent  breakdown".  In¬ 
stead,  the  "degeneration  to  turbulence"  would  be  due  to  the  development  of 
stronger  interactions  between  wave  modes.  As  the  waves  increase  in  amplitude 
and  slope,  the  interactions  between  modes  would  increase.  When  such  inter¬ 
actions  become  sufficiently  "promiscuous"  (to  use  his  apt  term)  one  would  then 
have  a  turbulent  cascade.  More  precisely,  Phillips  considered  wave  interactions 
that  involve  a  "resonant  triad".  The  latter  involves  a  three-wave  interaction  that 
satisfies  certain  special  conditions.  In  this  case,  he  showed  that  the  interaction 

time,  T.,  is  of  the  order 
1 

T.  -  (kt  o1k2  o2)‘l/2  (1) 

where  the  k's  and  ar's  refer  to  wave  numbers  and  fluid  particle  velocities,  respec¬ 
tively.  The  subindices  refer  to  two  of  the  three  particular  waves  in  the  resonant 
triad  (energy  goes  into  the  third  wave  from  waves  numbered  1  and  2).  As  the 
wave  amplitudes  increase,  so  also  do  the  particle  velocities;  hence,  from  Eq.  (1) 
T.  decreases.  This  decrease  of  T.  is  a  measure  of  the  increase  in  the  strength  of 
the  interaction.  The  shorter  T\  becomes,  the  less  wave -like  do  the  waves  in 
question  become  in  the  sense  that  they  do  not  obey  a  dispersion  relation.  In  other 
words,  if  a  mode  interaction  is  of  sufficient  strength  to  cause  T^  to  become  com¬ 
parable  to  a  wave  period  (inverse  H  runt -Vaisala  frequency),  the  wave  does  not 
propagate  but,  instead,  is  a  short-lived  local  entity  better  known  as  a  "turbulent 
eddy".  In  the  case  of  turbulent  eddies  there  is  no  longer  the  requirement  for  a 
resonant  triad  in  order  to  have  interaction.  The  following  are  two  of  the  crucial 
distinctions  between  waves  and  turbulent  eddies;  (a)  waves  obey  dispersion  re¬ 
lations;  eddies  or  turbulent  modes  do  not,  and  (b)  waves  last  a  long  time  and 
propagate  whereas  the  opposite  holds  for  eddies. 

Ideally,  a  wave  has  infinite  duration  (T.  =  x)  and  will  "linearly  superpose" 
with  other  waves  (that  is,  no  non-linear  interactions).  In  the  real  world  such 
requirements  can  only  be  approximated.  In  what  follows  we  will  not  use  so  ideal¬ 
ized  a  definition,  instead  we  shall  use  the  term  "waves"  in  a  sense  that  allows 
finite  duration  (such  that  the  decay  time  is  much  greater  than  a  buoyancy  period). 
We  will  also  use  "waves"  in  a  sense  that  allows  for  weak  non-linear  interactions 
such  that  T.  is  much  greater  than  a  buoyancy  period. 


Ideally,  an  <  .My  would  lose  its  identity  within  one  wave  oscillation  period  and 
within  one  v,  iv<  length  due  to  its  strong  interactions  with  other  modes  of  oscilla¬ 
tion.  This  is  presumed  to  be  a  good  description  of  turbulence  in  the  real  world. 
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>ns>d>- rations  can  shed  new  light  on  the  physical  meaning  of  the 
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t  .  (See  also  the  review  of  Phillips-".  )  This  buoyancy 
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(2) 


•.a*  ipatiou  rate,  the  buoyancy  frequency,  and  C  is  a  constant  of 

.  .  •.  deg  to  !  .umli  y,  this  length  separates  the  scales  of  the  buoy- 

t  i>  m  i  from  the  inertial  subrange,  the  latter  of  course 
: .  s  -tiles.  We  now  derive  ltq.  (2)  from  an  entirely  different 
i:  -  ph.v.-i  approach  than  given  by  I-umley. 

■ : . . »  :  i.  ■  ,  suppose  we  consider  the  ease  of  Kq.  (1)  where  k  =  k.,  =  k 

- 1  1  *- 

•  1 1  v.  \ls-j,  w  (  take  the  ease  where  T.  -  ,  which  is  to  say  that  we 

. lie  ■  .  li-i!  me  case  where  the  interaction  time  is  of  the  order  of  the  min- 
-aid.  This  is  the  case  >.  lu  re,  according  to  the  previous  section,  the 
.  ■■  :iuo[ ; ,  .•  a\  es  an  turbulent-like  in  their  pr  omiscuity  of  interaction. 

1  q.  ( 1 ),  w  e  have 

\  '3) 

i . 


>t‘u  ■  -if  the  wave  mode  in  question  gives  up  all  its  kinetic  energy  to  a  mode 

:  id  e  ra  ■  ‘  ion  cascade  ji.  on<  period  of  its  oscillation.  Now,  the  cascade  must  even- 

■  i : list  energy  :>t  lh<  molecular  ie'cci  at  tin  .  ;..a I : e;i!c  end  at  the  rate  e. 

•) 

I  - -king  v"  as  tla  1 : i i :c  t i •  energy  per  unit  of  mass,  and  assuming  the  cascade  to  be 
i  onservn live ,  we  arrive  at 


Phillips,  O.  M.  ( 1  ? ■ ' ; 7  )  On  tin  Holgiano  and  l.umley-Shur  Theories  of  the 
: .-M  .an-  \  subrange,  in  Atmospheric  t'url)>ilen<~i  and  Ilatlio  Wave  Prupaga- 
tioii,  A.M.  Miglom,  anil  V.  I,  Tatarsky,  Kds.  ,  NAUKA,  Moscow  ,  374  pp. 
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(4) 


Eliminating  v  in  Eq.  (3)  by  means  of  Eq.  (4)  and  solving  for  k,  we  arrive  at 


which  is  exactly  the  same,  in  effect,  as  Eq.  (2)  for  kg. 

This  argument  reveals  a  surprising  result.  The  buoyancy  length  does  not 
merely  divide  the  isotropic  from  the  anisotropic  turbulence.  In  addition,  it  also 
gives  the  minimum  wavelength  for  a  horizontally  propagating  gravity  wave.  Many 
new  questions  are  raised  by  the  above  considerations  about  the  nature  of  buoyancy 
turbulence.  For  this  reason,  these  issues  will  be  discussed  at  length  in  Section  8. 


I.  EXPERIMENTAL  CRITERIA  TO  DISTINGUISH  BETWEEN  WAVES 
AND  TURBULENCE 

2  0  ’’  1 

Busch  and  Stewart”  have  suggested  criteria  that  would  enable  one  to  deter¬ 
mine  experimentally  whether  or  not  a  given  field  of  fluctuations  is  due  to  turbu¬ 
lence  or  to  waves.  The  most  promising  of  these  criteria  is  based  on  the  fact 
mat  turbulence  causes  mixing  whereas  waves  do  not.  In  the  discussion  given 
previously,  the  strength  of  the  mode  interactions  was  the  criterion,  and  this 
raises  the  question  of  what  the  physical  connection  between  these  two  different 
criteria  might  be.  As  will  be  explained,  they  are  indeed  closely  related  physically. 
In  terms  of  practical  applications  concerning  pollution  transport  and  chemistry, 
as  well  as  optical  turbulence,  it  is  the  mixing  property  of  turbulence  that  is  the 
one  of  g i  eati-st  interest. 

Stewart”  1  suggested  that  one  simultaneously  measure  the  vertical  velocity 
flui  tuaiions  and  potential  temperature  fluctuations.  (In  place  of  temperature  one 
could  measure  any  other  scalar  quantity  such  as  the  concentration  or  mixing  ratio 
of  a  neutrally  buoyant  substance  provided  it  hail  a  significant  vertical  gradient). 
From  these  measurements  one  would  then  calculate  the  cohorcnct  ot  tween  these 


two  fluctuations  and  this  would  give  the  degree  of  vertical  mixing  taking  place. 

It  would,  in  fact,  give  the  vertical  transport  of  the  scalar  quantity  due  to  the 

turbulence.  In  the  case  of  waves,  this  flux  would  be  zero. 

2  0 

Busch  suggested  a  similar  idea,  but  he  introduced  a  useful  scale  depend¬ 
ence  as  follows.  Instead  of  coherence  he  suggested  that  one  calculate  a  cross  - 
spectrum  of  the  vertical  velocity  and  temperature  fluctuations.  This  spectrum, 
the  reader  will  recall,  has  a  real  part  (the  co-spectrum)  and  an  imaginary  part 
(the  quadrature  spectrum).  For  an  ideal  wave,  the  vertical  velocity  and  temper¬ 
ature  fluctuations  would  be  (±*/2)  out-of-phase  or  in  quadrature.  The  cross  - 
spectrum  would  be  purely  imaginary  and  there  would  be  no  net  vertical  transport. 
In  the  case  of  turbulence,  there  would  be  transport  or  mixing,  and  the  phase  of 
the  cross-spectrum  would  be  (±»),  that  is,  it  would  be  real  and  with  no  imaginary 
component.  The  phase  angle,  which  we  shall  call  |3,  will  in  general  depend  on  k, 
and  it  is  defined  as 

(3(k)  =  arc  tan  (  Re  $  ^(kf J  -  (6) 

where  $w,  i(k)  is  the  cross-spectrum  between  fluctuations  of  vertical  velocity 
w',  and  potential  temperature  O'.  In  general,  therefore,  for  those  values  of  k 
where  ji  =  90°  or  27  0°,  the  motion  would  be  wave-like;  and,  where  (J  -  0°,  or 
180°,  it  would  be  turbulent-like.  In  between  these  two  extreme  cases  the  motion 
would  be  neither  pure  turbulence  nor  pure  waves. 

Axford“^  tested  this  approach  experimentally  in  the  stratosphere.  He  found 
that,  for  the  cases  of  well-defined  trains  of  waves,  /i  =  90°  ±  10°.  Furthermore, 
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be  found  that  the  coherence  spectrum  lor  {[Im  ®wigi(k)l  +  (Re  ^w'g'^k)]  } 

normalized  by  the  individual  power  spectra  of  w'  and  o']  had  values  greater  than 
0.8.  This  second  finding  is  consistent  with  the  idea  that  a  wave-like  disturbance 
has  an  extended  periodic  pattern  in  space.  It  gives  a  second  test  for  waves  to 
use  in  conjunction  with  the  above  "(1  test".  In  the  cases  of  turbulent  motion  at 
the  scale  of  k,  Axford  found  that  -4r>°<  ;i(k)  <  +49°,  or  135°  <  p(k)  <  225°  were 
valuable  as  criteria. 

In  view  of  these  theoretical  and  experimental  findings  we  arrive  at  one  of 
the  conclusions  of  this  paper;  namely  it  would  be  desirable  to  apply  such  tests  to 
the  kind  of  data  described  by  Crooks  et  al,  ^  that  is,  those  upon  which  Figure  1 
is  based. 

23.  Axford,  lb  N.  (197  1)  Spectral  analysis  of  an  aircraft  observation  of  gravity 
waves,  Quart.  J.  Met.  Soc.  97:313-32  1. 


Other  suggestions  will  be  found  in  Busch.  We  only  mention  that  the 
individual  power  spectra  of  o'  and  of  w1  when  taken  in  combination  would  also  be 
helpful  in  distinguishing  between  turbulence  and  waves.  In  the  case  of  turbu¬ 
lence,  they  would  be  comparable  in  shape  and  magnitude.  In  the  case  of  waves 
they  would  differ  greatly  in  magnitude  with  the  o'  spectrum  being  much  smaller 
than  the  w’  spectrum,  especially  at  high  values  of  k. 


.-».  REMARKS  ON  THE  RELATION  BETWEEN  STRONG  MODE  INTERACTION 
\ND  “MIXING”,  AND  SPECTRAL  CONSIDERATIONS 


The  purpose  of  this  section  is  to  (a)  explain  the  physical  relationship  be¬ 
tween  strong  interactions  between  modes  of  oscillation  and  the  mixing  property 
of  turbulence.  In  the  process  we  shall  examine  the  physical  nature  of  turbu¬ 
lence,  going  well  beyond  the  previous  discussions;  and  (b)  examine  the  ambiguous 
nature  of  the  interpretation  of  the  power  spectrum  in  this  context.  The  latter 
seems  to  have  been  a  frequent  source  of  confusion  in  the  literature.  We  shall 
onciude  with  some  remarks  about  the  use  of  higher-order  spectra  for  fne  pur¬ 
pose  ;,f  dirt  .  tlv  measuring  strong  nteractions. 

A.  Ot  gu;  by  asking.  What  is  tne  connection  between  promiscuous  mode 
interaction  aim  mixing  (that  is  to  say.  mixing  in  physical  space)?"  Consider 
first  the  •  ast  of  idealized,  smaii -amplitude  waves  that  exhibit  perfect  super¬ 
position  or  non-mu  ruction  between  each  other.  1'ake,  for  example,  the  case  of 
a  fit  in  n  random  inti  rtuu  waves,  fan  these  cause  mixing?  A  single  wave  would 
not  causi  mixing  because  the  motion  would  oe  pel  iodic,  and  a  partieh  of  fluid 
w  ould  per uv iu  ally  return  to  a  fixed  location  in  space.  Two  w  aves  of  eomrm  usu¬ 
ry  tc  frequemy  would  obviously  have  the  same  effect,  but  the  time  between  returns 
•would  be  .01. g<  1  |i><  ir.g  :*  Uq  )  +  (p  m.,,))  where  the  w,  are  the  wave  frequencies 
and  p  aid  q  are  integers]  .  In  other  words,  the  motion  would  still  be  periodi.  . 

A  hat  .appt  r.s  in  the  ase  where  there  are  a  very  large  number  of  waves  w  ith 
,'c  ,  ; . -  male  ire  quem  les  '  .-nice  we  assume  omplete  supt  rpo.snioii,  any 
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is.  How  do  strong  inode -interactions  give 
If  a  wave  mode  decays  in  a  very  short  time, 
a  n  he  iei'l  in  their  initial  positions.  A  random  superposition  of 

ip.igating  modes  (that  is.  <  idles)  would  therefore  be  expected 
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Figure  2.  Vortex-Stretching  Mechanism  in  a  Turbulent  Cascade  (after 
Tennekes  and  Lumley^^).  The  stress  field  is  given  by  the  arrows. 
I.arger-scale  deformation  gives  energy  to  the  small-scale  vortices  by 
stretching  them.  This  causes  the  small-scale  vortices  to  increase 
their  spin  rotation  due  to  conservation  of  angular  momentum 
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the  cascade.  The  discussion  is  continued  in  an  article  by  Tennekes  where  he 
gives  a  convincing  argument  that  shows  that  a  large  eddy  breaks  up  into  two 
"daughters,  *  each  of  which  is  about  half  as  big  as  the  mother".  He  then  adds, 
"The  evolution  of  and  interaction  between  eddies  involve  non-linear  mixing,  both 
in  coordinate  space  and  in  wave  number  space.  "  As  was  already  mentioned,  the 
former  is  due  to  the  random  rates  of  strain  caused  by  the  cascade  process.  It 
should  now  be  clear  why  strong  mode-interaction  implies  mixing.  The  defor¬ 
mation,  since  it  separates  nearest  neighbors  of  fluid  particles,  is  the  cause  of 
the  mixing.  The  deformation,  in  turn,  is  due  to  mode  interactions. 

The  above  sort  of  cascade  picture  led  Kolmogorov  to  his  famous  -5/3  spec¬ 
trum  by  means  of  a  dimensional  argument  as  follows.  Since  the  power  density 

2 

spectrum  of  the  velocity  fluctuations,  ^(k),  has  the  dimensions  of  v  /k  or 
l  I./j  l  T*~|  ;  and  since  we  assume,  by  virtue  of  the  conservative  cascade,  that 
't(k)  depends  only  on  e  and  k  with  dimensions  [  l,  J  [  T  ]  and  L  respectively, 
it  follows  that 

-Mk)  k'/53  (7) 


That  is,  one  large  vortex  feeds  two  smaller  ones  until  the  large  one  is 
exhausted.  Onsager-^  envisioned  the  same  factor  of  2. 

Recently,  two  new  mathematical  approaches  have  appeared  that  shed  light 
on  some  of  the  classic  problems  of  this  subject.  The  interested  reader  may  con¬ 
sult  the  Appendix  for  a  brief  sketch  of  these. 

26.  Onsager,  L.  (1P4P)  Statistical  hydrodynamics,  Nuovo  Cimento  6(Ser.  P, 

Supp.  2  ):27  f.i-287. 


16 


where  a  is  a  constant  of  order  unity.  Note  that  in  this  derivation  it  was  not  neces¬ 
sary  to  assume  locality  of  interaction  in  k-space.  Rather,  all  that  is  used  is  the 
assumption  that  e  does  not  depend  on  k. 

We  now  turn  to  what  is  perhaps  the  key  to  the  confusion  that  leads  to  state¬ 
ments  like  "turbulence  is  a  field  of  random  waves";  and  "if  the  spectrum  has  -5/3 
slope,  this  means  ttiat  it  represents  cascade  or  inertial  range  turbulence".  The 

key  issue,  in  my  opinion,  lies  in  the  interpretation  of  the  power  spectrum. 

2  ;') 

I’ennekes  brings  this  out  very  clearly.  There  are,  in  fact,  two  possible  inter¬ 
pretations.  In  the  case  of  waves  one  can  imagine  the  Fourier  components  as 
representing  individual  waves  propagating  through  the  entire  sample  of  fluid  under 
consideration.  This  is  a  very  natural  wav  to  consider  Fourier  components  because 
this  is  the  way  the  subject  is  taught  in  courses.  Or:  tire  other  hand,  when  one 
considers  turbulent  motion,  one  is  not  looking  at  extended  waves  but,  instead,  at 
short-  'ived.  highly-localized  modes  of  disturbance  or  eddies.  One  is  considering 
a  large  o. lection  or  assemblage  of  these  entities  (approximating  an  ensemble) 
ard  the  Fouriei  <  omponent  does  not  refer  to  a  single  entity  su<  h  as  a  wave. 

Instead,  as  Tennekcs  puts  it,  "\Ve  conclude  that  a  Fourier  coefficient  should  not 
oe  contused  with  an  eddy.  Fourier  coefficients  are  associated  with  an  ensemble 
:>f  eddies,  and  when  one  ana  vzes  the  dynamics  of  turbulence  in  wave  number 
sp.ua  one  sno  'c.  use  tire  plural  otuv:  (Ui<  oueettve  of)  eddies  at  or  near'  a  cer¬ 
ium  wave  number  are  mvoiveu.  Flit-  reader  stiouici  consult  his  work  lor'  his  use 
>1  wave  paeket  and  quantum  mechanical  analogies,  which  are  very  illuminating. 

*  »t  <  o'.rsi  it.  vte  a  ot  Hie  preceding  discussion,  there  is  a  complete  range  ol 
lynamtes  note,  v  avis  to,  slum  we  saw,  lurb-undulenco",  and  finally  to  turbulence; 
i:i  :  a  pow  <  ;  spei  trurn  alone  says  nomine  about  inis  aspect  of  the  situation.  As 
was  already  explained,  a  .  ross-spew  trum  can  bring  out  the  difference  between 
waves  am,  timouiem  i  .  I  ius  mi  too  c pends  an  o-tn  upon  the  mixing  property 
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frequencies  will  be  found  in  Dewan.  Such  spectra  measure  non-linear  inter¬ 
actions  between  "wave-modes".  In  the  absence  of  such  non-linear  interactions 
the  bi-spectrum,  for  example,  would  be  zero. 

Unfortunately,  the  idea  of  using  higher-order  spectra  to  measure  the  transfer 

rate  of  energy  between  modes  of  oceanic  internal  waves  has  problems  associated 
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with  it  that  render  it  impractical  at  this  time  (as  explained  by  McComas  ).  On 
the  other  hand,  the  application  of  higher-order  spectra  to  atmospheric  turbulence 
has  yet  to  be  explored.  In  the  case  of  turbulence,  there  is  every  reason  to  be 
optimistic  about  the  usefulness  of  higher-order  spectral  measurements.  For  the 
present,  then,  only  mixing  in  coordinate  space  is  available  for  measuring  "degree 
of  turbulent  interactions"  in  a  fluctuation  field.  However,  higher-order  spectra 
may  very  likely  play  a  role  in  the  future. 


6.  TWO  DIMENSIONAL  TURBULENCE 

The  term  "two-dimensional  turbulence"  is  well  known,  and  was  first  discussed 
by  Onsager36  in  terms  of  a  "reverse  cascade".  More  recently,  Kraichnan3^  dis¬ 
cussed  this  two-dimensional,  reversed  cascade  (large  k  to  small  k)  of  energy, 
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and  he  showed  that  it  would  have  a  k  ‘  '  dependence.  Leith  and  Kraichnan 

subsequently  showed  that  in  two-dimensional  turbulence  there  is  also  the  possi- 

-3 

bility  of  a  forward  cascade  of  enstrophy  that  would  have  a  k  dependence  for  the 
spectrum  (see  also  Charney33).  Both  Stewart31  and  Gage34  have  proposed  that 
some  of  the  observations  of  k  '’^3  one-dimensional  spectra  in  the  atmosphere  may 
be  the  result  of  a  reverse  two-dimensional  cascade.  This  leads  directly  to  the 
question  of  whether  or  not  Figure  1  could  be  explained  along  such  lines. 

An  experimental  approach  to  the  answer  to  this  question  would  be  to  invent  a 
modification  of  the  Busch  and  Stewart  approach  that  would  render  it  appropriate 
to  two-dimensional  turbulence.  We  now  attempt  to  do  this.  In  the  previous  case, 
we  had  three-dimensional  turbulence  in  which  was  imbedded  a  scalar  quantity,  0 
(the  potential  temperature),  which  had  a  vertical  gradient.  Vertical  transport  due 

20.  Dewan,  K.  M,  (lOMP)  Nonlinear  Cross -spectral  Analysis  and  Pattern 
Kecognition,  A F (  K L-O  0  -  002  fi ,  AD  fi 87 085 1 . 

30.  McC  omas,  (  .11.  (1078)  Private  communication. 

31.  Kraichnan,  H.H.  (1007)  Inertial  ranges  in  two  dimensional  turbulence, 

Phys.  of  Fluids  10:1417-1423. 

32.  Leith,  ('.  K.  ,  and  Kraichnan,  H.H.  (1072 )  Predictability  of  turbulent  flows, 

■L  Atmos.  Sci.  20:1041-1088. 

33.  ('barney,  J.G.  (107  1)  Geostrophic  turbulence,  J.  Atmos.  Sci.  2  8:1 087  -  1 005 . 

34.  Gage,  K.S.  (107  0)  Kvidence  for  a  k  law  inertial  range  in  mesoscale  two- 

dimensional  turbulence.  ,1.  Atmos.  Sci.  3(i:105o-l(i54. 


to  mixing  provided  the  criterion  for  the  presence  of  turbulence.  In  the  case  of 
(horizontal)  two-dimensional  turbulence,  the  analogous  situation  would  have  to 
prevail,  that  is,  there  would  have  to  be  a  scalar  quantity  6  present  that  had  a 
horizontal  gradient.  Two-dimensional  turbulence  would  then  be  signaled  by  trans¬ 
port  via  the  horizontal  velocity  fluctuations  in  the  direction  of  V@,  and  one  would 
merely  substitute  that  velocity  component  for  what  we  called  w  previously.  The 
test  would  be  formally  identical  to  what  was  given  previously.  Whether  such  a 
test  would  be  practical  remains  to  be  seen;  for,  in  the  previous  situation  vertical 
mean  motion  was  negligible.  In  contrast,  there  would  be  horizontal  mean  motions 
that  would  have  to  be  eliminated  from  the  analysis. 

Theoretically,  there  are  two  problems  with  a  two-dimensional  turbulence 
explanation  of  Figure  1.  The  first  is  the  question  of  how  energy  could  be  supplied 
to  the  small-scale  end  of  the  spectrum.  The  second  is  the  question  of  how  to 
explain  the  observation  that  the  vertical  velocity  fluctuations  have  the  same  mag¬ 
nitude  and  power  spectrum  slope  as  the  horizontal  velocity  components  in 
(Figure  l). 

There  is  one  complication  in  this  experimental  criterion  for  turbulence  that 
was  not  pointed  out  by  Busch,  but  that  we  should  mention  for  completeness.  The 
cross-spectra  involved  are  based  on  one-dimensional  trajectories,  that  is,  we 
deal  with  one-dimensional  spectra.  It  is  well  known  that  such  spectra  are  sub¬ 
ject  to  spatial  aliasing  effects.  This,  in  turn,  could  introduce  an  ambiguity  into 
the  meaning  of  -;(k),  or  phase,  if  more  than  one  wave  gave  the  same  k  component. 

In  principle  this  could  be  remedied  by  making  use  of  higher  dimensional  spectra; 
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however,  the  notable  success  of  Axford  with  this  technique  suggests  that  my 
concern  is  premature. 

7  Till  Vi  \\  V  \SC\I)E  explanation  FOR  A  k  5/:l  SPEOTR l  M 

Mct'omas  1  anti  Mct'omas  and  Bretherton  have  investigated  internal  wave 
interactions  and  the  transfer  of  action  in  k-spaee.  Their  work  is  primarily 


33.  Mct'omas,  t'.H.  (1U77)  Equilibrium  mechanisms  within  the  oceanic  internal 
wave  field,  J.  Phys.  Oeeanog.  7:838-845. 

3*i.  Mct'omas,  C.M.,  and  Bretherton,  F.  P.  ( 1  r*7 7 )  Resonant  interaction  of 
oceanic  waves,  J.  Geophys.  lies.  82  : 1 3 f >7  —  1 4  12. 

-  ~~~  J  u 

Wave  action  is  defined  as  wave  energy  divided  by  frequency  (see  Phillips  ). 
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built  upon  the  former  researches  of  Phillips  and  Hasselman  *  among  others. 

Briefly,  their  treatment  assumes  that  the  waves  are  essentially  linear  and  that 

1  8 

the  interaction  between  modes  is  weak.  Phillips  showed  that  there  is  a  way  by 
which  internal  waves  can  interact  in  a  resonant  manner  provided  that  the  temporal 
frequencies  and  wave  numbers  of  a  triad  of  waves  obey  a  certain  condition.  This 
condition  is 

kl  1  k2  =  k3 

u.  ±  ujg  =  u3  (8) 

(In  the  case  of  surface  waves  there  must  be  a  four-wave  interaction.  ) 

The  physical  mechanism  (by  which  energy  can  be  transferred  from  two  inter¬ 
acting  modes  into  a  third  mode)  resembles  forced  simple  harmonic  oscillation. 
Two  modes  interact  by  means  of  the  convective,  non-linear  term  in  the  equations 
of  motion  and,  as  a  result,  they  generate  sum  and  difference  frequencies  [see 
Icq.  (8)|  .  The  third  wave  (number  3),  is  assumed  to  have  both  a  spatial  and  tem¬ 
poral  resonance  with  this  "beat  frequency  phenomenon"  caused  by  the  interaction 
between  waves  numbered  1  and  2.  Wave  three  is  thus  "driven"  in  a  manner 
analogous  to  a  simple  harmonic  oscillation  with  a  forcing  term  on  the  right-hand 

side  of  the  equation. 
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McComas  '  shows  that  among  all  the  possibilities  for  Eq.  (8)  to  hold,  there 

are  three  distinct  limiting  configurations  compatible  with  the  dispersion  relations; 

and  these  he  labels  "induced  diffusion",  "elastic  scattering",  and  "subharmonic 

instability".  These  allow  a  simplified  and  physical  way  to  further  understand  the 

transfer  between  wave  modes.  He,  for  the  purpose  of  this  paper,  provides  us 

with  a  proof  that  a  wave-wave  energy  cascade  can  exist.  Beyond  this  one  fact 

our  present  work  needs  no  further  details. 

Parenthetically,  it  should  be  mentioned  at  this  point  that  the  thrust  of  the 

3  G 

work  by  McComas  and  Bretherton'  ’  is  primarily  to  provide  a  theory  for  the 
"tiarrett-Munk"  Spectrum,"  '  which  represents  a  universal  wave  spectrum  to 
be  found  associated  with  internal  waves  in  the  ocean.  A  desirable  achievement 

37.  Masselmann,  K.  (i!>r>2)  On  the  non-linear  energy  transfer  in  a  gravity  wave 
spectrum.  Part  1,  .K  fluid  Mech.  12 :48 1-500. 

33.  Masselmann,  Is.  (1013)  On  the  non-linear  energy  transfer  in  a  gravity  wave 
spectrum.  Part  2  and  Part  3,  J.  Fluid  iVlcch.  15:273-281,  385-308. 

3'i.  Garrett,  (  .,  and  Munk,  W.  (1072)  Space-time  scales  of  internal  waves, 
Cieophvs.  Fluid  Dyn.  2:22  5-284. 

to.  Oarrett,  ,  and  Munk,  W.  (1075)  Space-time  scales  of  internal  waves:  a 
progress  report,  J.  Geophys.  Res.  80:291-207. 


would  be  to  explain  the  atmospheric  data  such  as  that  exemplified  by  Figure  1. 

Far  from  being  "universal",  such  data  are  found  only  on  rare  occasions  in  which 
the  wave  intensities  are  so  great  that  they  were  originally  identified  as  turbulence. 
Such  disturbances  are  found  on  the  order  of  5  percent  of  the  time.  Furthermore, 
the  disturbance  is  found  to  be  confined  to  altitude  thicknesses  of  1  to  2  km  (as 
shown  by  Crooks  et  al*'  lfi).  At  present  it  is  not  clear  if  our  conclusions  have 
any  application  to  the  Garrett-Munk  (GM)  spectra;  but  this  cannot  be  ruled  out  at 
this  time.  On  the  other  hand,  they  do  appear  to  have  relevance  to  Figure  1. 

In  Dewan,  an  alternative  existence  proof  of  the  possibility  of  an  energy 
wave  cascade  is  given.  There,  the  interaction  between  wave  scales  is  explained 
not  by  resonance,  but  by  the  effect  of  deformations  caused  by  large  waves,  which 
in  turn  can  "feed"  small  waves.  The  next  step,  if  one  were  to  work  in  parallel 
with  the  bizarre  history  of  turbulence  theory  (that  is,  Richardson’s  famous  poem 
followed  by  Kolmogorov's  dimensional  argument)  would  be  to  create  a  new  poem, 
"Big  waves  have  little  waves  that  feed  on  deformation  and  little  waves  have  lesser 
waves,  to  turbulent  dissipation".  Needless  to  say,  the  poem  is  not  a  prerequisite 
to  making  a  dimensional  argument  nor  does  it  reflect  the  resonant  mode  inter¬ 
action  approach  of  this  paper.  Instead,  it  reflects  an  alternative  approach  given 
in  Dewan.  ^ 

In  one  sense  there  is  a  very  close  parallel  between  ne  wave  cascade  theory 
17  35 

given  in  Dewan  and  that  given  in  McComas.  The  latter  involves  interactions 

with  time  scales  given  by  Kq.  (1)  or  its  simplification  given  by  T^  —  (kv)  \  which 

was  used  in  Section  3.  In  contrast,  the  treatment  given  in  Dewan17  involves 

production  terms  of  the  form  uw  9U/8z  where  U  represents  average  larger-scale 

velocity  fluctuations,  and  u  and  w  the  smaller  scale  ones.  These  production 

terms  represent  transfer  from  larger  to  smaller  scales.  In  this  case,  however, 

T.  =  (9U/3z)  .  However,  in  a  wave  cascade,  U  is  oscillatory  and  its  shear  is 

1  —  «  1 

given  by  Uk  thus  (T.  )  is  physically  the  same  quantity  as  that  given  in  resonance 
1  o 

theory.  Instead  of  e  =  v"  (vk)  for  resonant  mode  interactions  (where  t  is  the 

J  W  'V 

transfer  rate)  one  would  have  f  =  v“  dv/dz  which,  in  effect,  is  the  same  quantity. 

The  cascade  in  question  must  involve  transfer  from  the  large  to  the  small  scales. 

This  is  well  known  because  turbulence  itself  (at  the  small  scale)  is  dissipative 

and  ultimately  converts  motion  into  heat.  It  cannot  represent  the  original  source. 

The  latter  must  be  at  large  scales.  True,  turbulence  can  generate  gravity  waves; 

however,  it  itself,  must  be  generated  by  shears  that  arise  from  larger  scale 

17 

fluctuations.  For  more  discussion  of  this  the  reader  should  consult  Dewan. 

In  any  case  we  now  make  the  assertion  that  the  wave  power  spectrum  depends 
only  on  k  and  (which  is  the  dissipation  rate  for  the  waves  due  to  their  inter¬ 

actions).  in  order  to  calculate  the  actual  kinetic  energy  dissipation  rate  per  unit 
/olume  one  must  multiply  by  ij ,  the  mass  density. 


We  are  now  in  a  position  to  write  down  the  expression  for  the  wave  cascade 

one -dimension  spectrum.  One  need  only  repeat  the  dimensional  argument  of 
9 

Kolmogorov'  given  previously  to  arrive  at 

*  _  2/3  .  -5/3  ... 

=  °W  eW  ^  •  (9) 

Here  the  subscript  W  signifies  waves. 

One  could  object  that,  unlike  turbulence,  wives  must  obey  a  dispersion 
relation  that,  in  turn,  must  depend  on  the  buoyancy  frequency.  In  its  simplest 
form,  this  dispersion  relation  is 

u  =  Ng  cos  e  *  (10) 

4  1 

where  f)  is  the  angle  (q  £  1/2 it)  between  the  horizontal  direction  and  k  (Turner  ). 
Thus,  how  can  one  assume  that  $(k)  does  not  depend  upon  Ng?  My  answer  is 
that  we  presently  are  assuming  a  cascade  in  steady-state  and  in  which  negligible 
mixing  is  taking  place.  (We,  like  McComas,  assume  that  the  interactions  are 
weak  in  comparison  to  turbulence.  )  That  is  to  say,  tyj  output  would  equal  the 
input  so  long  as  energy  is  not  given  up  in  the  form  of  the  potential  energy  due  to 
mixing  in  a  stratified  environment.  This  implies  that  buoyancy  does  not  affect 
f^..  It  is  therefore  reasonable  to  suppose  that  a  change  in  Ng  would  affect  only 
the  amplitude  of  the  waves  in  such  a  way  that  the  kinetic  energy  for  a  scale,  k, 
would  remain  unaltered.  In  other  words,  our  dimensional  argument  is  physically 
permitted. 

The  other  facts  that  make  the  hypothesis  of  a  wave  cascade  a  useful  one  are 

(a)  the  experimental  values  of  the  spectra  in  fact  give  us  a  k  1  '  ,  dependence, 

(b)  as  will  be  seen,  can  be  related  to  observation,  and  (c)  one  can  experi¬ 
mentally  test  the  "wave  field"  hypothesis  via  the  previously  described  "Busch 
Test"  to  verify  that  it  is  indeed  a  wave  field.  In  other  words,  this  approach  is 
testable. 

Finally,  we  now  consider  the  theoretical  shape  of  a  spectrum  that  extends 
through  ranges  of  k  including  both  the  waves  and  turbulence.  This  requires  that 
we  obtain  the  answers  to  the  following  three  questions:  (a)  "At  what  value  of  k 
is  the  transition  between  'f’^fk)  and  ^(k)  (the  turbulence  spectrum)?,  (b)  "in 


■■  9  2  2  2 

More  generally,  w“  =  N.,  cos  0  +  w.  sin  n  where  *j.  is  the  inertial  frequency 
(see  Garrett  and  Munk-^»  40).  1  1 

41.  Turner,  J.  S,  (1P73)  Buoyancy  Effects  in  Fluids,  Cambridge  University 
Press,  307  pp. 
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what  way  is  ew  related  to  eT  (turbulent  dissipation)"?  and,  (c)  "What  are  the 
numerical  values  for  a ^  and  ?" 

The  first  question  is  easily  answered  on  the  basis  of  the  previous  discussion. 
The  "buoyancy  length"  derived  above  gives  us  the  boundary  between  horizontally 
propagating  waves  and  isotropic  turbulence.  While  this  "boundary"  is  a  transition, 
and  is  hence  smooth,  it  is  given  by  kB  where 


(see  Eq.  (5)]  where  c  is  a  constant  of  order  unity. 

The  answer  to  the  second  question  (regarding  the  relation  between  the  wave 
and  turbulence  dissipation  rates  e^,,  and  trp  respectively)  is  obtained  from  the 
assumption  that  the  cascade  is  conservative.  Energy  conservation  gives 

p  ~  P  eT^T  (12) 

where  V^,  represents  the  volume  occupied  by  the  waves  and  V^.  that  of  the  turbu¬ 
lence.  The  reason  V^,  and  VW  are  not  equal  in  a  stratified  medium  is  that  turbu¬ 
lence  in  the  latter  is  known  to  occur  in  thin  layers  occupying  a  small  fraction  of 
the  total  volume.  Thus  we  arrive  at 


(13) 


Since  V^,  >  V.^,,  it  follows  that  ew  <  e  .  Knowing  ew,  one  could  calculate  if 
the  volume  ratio  were  known.  It  should  be  noted  that  if  one  were  to  ascertain 
'f.j,  by  means  of  the  spectrum  (assuming  at  that  point  that  and  were 
known)  then  one  would  have  a  measure  of  V^,/V^.  that  would  in  turn  have  some 
pragmatic  utility  in  problems  concerning  pollution  transport  and  optical  trans¬ 
mission  through  turbulence.  Specifically,  V -p/V\y  would  be  a  measure  of  the 
volume  occupied  hv  turbulence  relative  to  the  overall  volume  of  the  disturbed 
region. 

rile  third  question  (concerning  the  universal  constants  of  order  unity,  and 
<>,.  )  can  only  be  answered  by  experiment.  One  already  knows  that  a™  is  around 
1.  '>  (Tennekes  and  l,umlcy“  )  but  it  depends  somewhat  on  the  velocity  component 
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in  question.  (One  does  not  yet  know  the  value  of  but,  according  to  Bond, 

there  is  a  PC  percent  a  priori  chance  that  it  lies  between  0.  l  and  10.  ) 

Figure  3  shows  a  schematic  representation  of  this  spectrum  assuming  that  the 

measurement  is  made  from  a  probe  within  a  turbulent  layer. In  this  figure  no 

account  has  been  taken  of  the  spatial  aliasing  effects  to  be  expected  in  the  case  of 

4  3 

one -dimensional  spectra.  Gifford  showed  ct ii  CXSiiJ  pie  where  there  was  a  —5/3 
slope  that  continued  to  wavelengths  three  to  five  times  larger  than  would  appear  in 
a  spectrum  without  aliasing. 

To  il lustra t*.  this  effect,  let  us  consider  a  specific  numerical  example.  Sup¬ 
pose  that  kB 1  =  aO  m  (a  typical  value  for  HU  AT  spectra  such  as  in  Figure  l). 

In  this  case,  A.  =  3  1 3  m.  If  aliasing  were  to  extend  the  —  f; / 3  slope  to  scales  live 
times  larger  than  A  that  would  take  it  out  to  1,  >7 . ■  m.  f  igure  1,  however,  has 
the  slope  of  -”>/3  extending  to  about  one  order  of  magnitude  beyond  this  length. 

It  would  be  very  useful  il  velocity  spectra  and  cross  spectra  [l.q.  (r;)J  wire 
measured  in  stable  fluids  in  order  to  test  the  idea  that,  at  wavelengths  much 
iarger  than  waves  exist,  and  much  less  man  A.,,  turbulence  exists,  and  that 
both  can  have  a  -  3  si  opt  but  with  t?  ^  *  >  p  (.-.e<  i  q.  ( 1  :>  >  |  .  l.quation  (('•),  which 
is  wiiat  l  b.avi  cairn  :  nu  Bus- -it  test,  .  i  pi  ui.-.i  rn  between  w aves 

and  turimience.  ip  t  tins  pom;,  wt  na  .  ■  lett  out  hue. an  umirdeiice. 

a  on  mi  nut  mol  mm  in  oi  tv  i  o  if  if  \N/,i  <n  n  kiu  u  \n 
V  N 1 1  \  M  \  I  t  I  I  N(. Mi  '  toil  Kit  >  N 

i n  tin 1  i ■  i •  >  v  ;•  '■ :  .  i  •  ■'  .  :  :  •  :  •  a,.  a •  oi  wavt  s  a  ill': 

iua>:  ,  . .  :  :  ;  i.-, ;  ts  I ;  i  >: 

tile  hori/onta .  It;-  ;  ■  •  ■  ■  .  >i  u  ■  t  •.  hori.'oi  - 

tailvl,  tin  >  i  •  •  ;t.  . 

Ac  ‘.;U’  ;;  i  ■  '  '  •  >:  ■  r .  .  i  ■.  <  :  ; 

lu  t  •.«  \  \  ■  -  . •  • ;  •  i  ,  •  .  .  .  •  .  .  >  .  :  .  \  ■  :  . ■  •  : . i  • 
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Jsing  the  fact  that  (v/i  )  is  of  the  same  order  as  the  mean  vertical  shear 
(Reference  24),  we  can  write  the  Richardson  number,  in  the  form 


R  =  - q 

1  (v'f  r 

o 


A  typical  value  for  R^  in  cases  where  shear  induced  turbulence  exists  is  (l/4). 

Using  this  it  is  easy  to  obtain 

i.-'1  =  '  (2.88)  .  (20) 

l-j  o 

Hms,  within  a  constant  of  >rd(  r  unity,  the  outer  length  and  buoyancy  length  are 

i  nun,’. 

This  raises  some  interesting  questions.  Usually  1  defines  the  scale  at 
v ! 1 1 1  h  the  turbulence  is  generated.  Hut  k.,1  is  the  minimum  scale  for  buoyancy 

I  j  - - - 

subrange  turbulence.  How  ttien  -an  there  be  buoyancy  turbulence  if  the  energy 
■ssi  :eies  to  smaller  scales'’  1  ould  it  be  the  case  that  buoyancy  turbulence  must 
be  generated  by  a  larger  scale  horizontal  shear  ?  This  is  not  an  unreasonable 
i.oKsibiiity  since,  as  w  ,  ha1.'  seen,  this  kind  of  turbulence  involves  horizontal 
eddii  s  (espei  laliv  at  large  vnint  s  of  a).  Another  possibility  that  comes  to  mind 
is  that  a  reverse'  cascade  may  lit-  involved  in  buoyancy  turbulence.  It.  is  well 
Une.wti  that  two-dimensional  turbulence  can  involve  a  reverse  cascade,  and  for 
large  i  ,ve  hat  e  'W  .n  that  buoyancy  turbulence  becomes  two  dimensional. 

;  .i  aging  tlic.se  questions  for  future  research,  we  now  turn  to  the  use  of  Kq.  (In) 
■is  ,,  .  to  distinguish  turbuh  nt  from  wave  fluctuations.  It  has  the  unfortunate 

i r ; i  v1'.-,  k  that,  it  requires  knowledge  of  the  value  of  0  for  the  k  under  consideration, 
b.  ;■(•!.  ip. i  -u  •!  information  .von!  ;  or  ;iv:i ilabu.  if  the  complete  three -dimensional 
■  ’  ■  av<  -nur.bf  r  c  •  trim  wer*  ..-a  -mr'-d.  it  hiie  su>  h  ine-isur'-ments  might 

.  i : i  m  i  .n  ornlorv,  •(;.  v  -in  ;sua!lv  impossible  to  ordain  in  atmospheric 

.  1  r  *(,.  . , « i  i  •  •  i  ■  •  i;  1 1 1  .u  r  ills.-  nssior.  suggests  that  tin  rt  may  be  signal 

n<  pi  b,  *;m.  ,f  da  :  ■  n  'nations  can  ne  used  as  a  wav  turbuii  nee 


>.:•,'/<  'n  i  :’i.  dug  situation  that  art  trapy i.avt  l*  -  ili 

•  .  s,  a  n  :>  du  i.  leaving  ut  him  only  tin  nisotropu 


turbulence  of  large  X  and  the  horizontally  propagating  wa/es.  As  we  have  seen, 
such  turbulence  would  involve  mainly  horizontal  motion.  Thus  vertical  motion 
would  be  due  to  waves  and  horizontal  motion  would  be  due  to  eddies,  making  the 
distinction  merely  a  matter  of  measuring  the  direction  of  the  velocity  fluctuations. 


9.  CONCLUSIONS 

Waves  and  turbulence  are  idealizations  for  two  extremely  different  types  of 
fluctuations  in  a  stable  fluid.  Waves  represent  the  extreme  of  linear  (super- 
posable)  dissipationless  motion.  Turbulence  represents  very  strong  mode-inter- 
actions,  mixing,  and  dissipation.  Waves  have  non-local,  long-lived  propagating 
patterns,  and  turbulence  has  local,  short-lived,  non-propagating  patterns. 

E  have  used  "waves"  in  this  paper  in  a  manner  that  departs  somewhat  from 
the  ideal  type  in  that  some  weak  non-linear  interactions  are  permitted.  In  any 
case,  the  cross -spectral  test  of  Busch  was  discussed  and  it  should  prove  useful 
in  distinguishing  between  waves  and  turbulence.  There  is  reason  to  hope  that 
bi-spectra  could  also  play  a  future  role  in  making  this  experimental  distinction. 
More  specifically,  the  "bi-coherence"  defined  in  MeComas30  could  be  used  to 
estimate  the  degree  of  turbulent  non-linear  interaction. 

Figure  1  presents  us  with  an  interesting  scientific  mystery.  On  one  hand, 
it  is  clear  that  the  inertial  turbulence  interpretation  is  ruled  out.  A  wave  cascade 
possibility  has  been  proposed.  Buoyancy  turbulence  or  two-dimensional  turbu¬ 
lence  (they  seem  to  be  related)  may  be  involved.  The  answers  to  the  questions 
raised  will  probably  necessitate  the  use  of  some  of  the  wave  turbulence  criteria 
proposed  in  this  report. 

Hopefully  this  work  will  prove  useful  in  the  practical  problems  of  pollution 

transport  and  optical  turbulence.  If  it  turns  out  that  the  cascade  hypothesis  is 

related  to  universal  wave  spectra  (Garrett  and  Munk3'1’ 4(*),  this  approach  could 

lend  to  better  ways  to  predict  shear  structure  giving  rise  to  turbulence.  If  it 

could  lend  to  a  better  understanding  of  the  intense  trapped  gravity  wave  structure 

of  the  type  presumably  observed  by  (  'rooks  et  al,  '  ’  ' fl  perhaps  this  could  lead  to 
•> 

ways  to  predict  high  conditions. 


Otic  mechanism  that  might  account  for  these  high-intensity  trapped 
could  involve  over-reflection,  for  example. 


■.Hi 
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Appendix  A 

New  Developments  in  Turbulence  Theory 


New  work  in  the  mathematical  approach  to  turbulence  enhances  our  under¬ 
standing  of  the  chaotic  nature  of  turbulent  mixing.  This  appendix  is  intended  to 
bring  this  new  information  to  the  reader's  attention.  The  major  unresolved  ques¬ 
tion  in  turbulence  theory  is,  "How  do  the  random  or  chaotic  motions,  which  sup¬ 
posedly  are  the  solutions  of  the  purely  deterministic  equations  of  fluid  motion, 
come  into  existence?"  The  new  findings  demonstrate  a  physical  mechanism  that 
makes  such  a  phenomenon  possible. 

It  is  very  well  known  that,  as  the  Reynolds  number  is  increased,  a  critical 
value  can  be  reached  where  the  laminar  flow  makes  a  transition  to  turbulent  flow. 
I'he  transition,  in  some  cases,  is  spread  over  a  range  of  Reynolds  numbers,  that 
is,  there  is  a  series  of  transitions  or  bifurcations  in  the  complexity  of  the  motion. 
As  Reynolds  number  increases,  complexity  increases  until  a  point  is  reached 
where,  suddenly,  the  motion  becomes  totally  chaotic. 

I’he  "breakthrough"  con-  ept  that  explains  the  chaotic  behavior  is  known  as  a 
"strange  attractor"  and  it  was  introduced  by  Lorenz,  '  *  in  the  context  of  the  Menard 
instability.  My  Fourier  decomposing  certain  variables  of  the  partial  differential 
equations  for  the  problem  he  obtained  an  infinite  set  of  ordinary  differential  equa¬ 
tions.  These  he  truncated  to  three  equations,  and  hence  his  work  represents  a 
severe  mutilation  of  the  actual  i  quations.  In  any  case,  it  led  to  the  first  discovery 

A  I.  I.orenz,  K.  N.  ( 1  f  *  ( ;  X )  Deterministic  nonperiodic  flow,  -I.  Atm.  Sei.  _M>:  i;iu- if  1 . 


A2 

of  strange  attractors.  Subsequently,  Ruelle  and  Takens  arrived  at  strange 

attractors  for  turbulence  theory  from  an  entirely  different  direction,  namely 

from  the  qualitative  theory  of  differential  equations.  They  did  much  to  explain 

both  chaotic  motion  and  also  the  Reynolds  number  bifurcation. sequence  described 

A3 

previously.  Most  recently,  the  work  of  Feigenbaum  provided  an  especially 

clear  way  to  understand  both  the  transition  to  increasing  degrees  of  complexity  of 

motion  as  well  as  the  subsequent  transition  to  total  chaos. 

In  order  to  understand  qualitatively  the  meaning  of  the  "strange  attractor" 

concept,  one  can  start  with  the  simpler  entity  known  as  an  attractor.  The  simplest 

examples  of  an  attractor  occur  in  works  on  non-linear  mechanics  (for  example, 

A4  A  5 

Andronov  et  al  or  Minorsky  ).  Consider,  for  example,  the  motion  of  an 
electronic  oscillator  as  modeled  by  the  van  der  Pol  equation.  Such  motion  is  con¬ 
veniently  examined  in  the  x-v  (that  is,  displacement  and  velocity)  or  "phase  plane". 
For  example,  if  the  oscillator  were  to  start  with  initial  conditions  near  x  =  0  and 
v  =  0,  the  trajectory  of  the  solution  would  form  an  outward-going  spiral  that 
asymptotically  approaches  an  ellipse  that  is  known  as  "limit  cycle".  The  latter 
trajectory  represents  an  equilibrium  between  input  and  output  of  energy  over  the 
cycle.  If  the  cycle  were  to  start  with  large  values  of  x  and  v,  the  trajectory 
would  spiral  inward  towards  the  limit  cycle.  The  limit  cycle  thus  appears  to 
attract  all  trajectories  no  matter  what  their  initial  conditions  might  be,  and  hence 
the  term  "attractor"  is  used. 

The  adjective  "strange",  when  employed  in  the  present  context,  designates  a 
key  aspect  of  the  explanation  of  chaotic  motion.  It  refers  to  the  fact  that  such  an 
attractor  is  of  fractional  dimension  or  is  a  "fractal"  as  is  described  at  great 

A 

length  in  the  book  by  Mandelbrot.'  As  he  explains,  a  fractal  curve  (or  at  least 
one  major  type  of  such  curve)  is  one  of  which  each  piece  is  a  reduced  scale 
version  of  the  whole.  Such  curves  have  the  property  that,  as  their  length  is 
measured  at  ever  higher  resolutions,  it  increases  indefinitely.  The  dimension, 

D,  is  given  by 

D  =  log  N./log  (l/r)  (Al) 


A2.  Ruelle,  D.  ,  and  Takens,  F.  ( 1  r>7  1 )  On  the  nature  of  turbulence,  C'ommun. 

Math.  Phys.  20: 167 - 1  f*2 . 

A3.  Feigenbaum,  M.J.  (1980)  Universal  behavior  in  nonlinear  systems, 

Los  Alamos  Science  Summer :4. 

A4.  Andronov,  A.  A.,  Vitt,  A.  A.,  and  Khaiken,  S.  F.  (190(1)  Theory  of  Oscillators, 
Pergamon  Press,  New  York,  815  pp. 

Af>.  Minorsky,  N.  (1947)  Introduction  to  Non-Linear  Mechanics,  Ldwards,  404  pp. 

A(i.  Mandelbrot,  B.  13.  (1977)  Fractals,  Form,  Chance,  and  Dimension, 

Freeman  and  Co.  ,  365  pp.  —  • 
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where  N  is  the  number  of  parts  in  a  whole  piece  of  the  curve  and  r  is  the  so-called 

similarity  ratio.  For  example,  consider  the  case  of  a  square  which  is  subdivided 

into  N  squares;  the  "similarity  ratio"  would  be  r  =  l/N1'  .  That  is  to  say,  each 

side  of  the  small  squares  would  have  this  ratio  to  the  length  of  the  original  square. 

Here,  one  can  obtain  D  =  2  from  Eq.  (Al).  But  D  =  2  dimensions  is  the  ordinary 

situation  for  plane  figures.  Thus  Eq.  (Al)  gives  the  usual  results  for  non-fractal 

objects.  The  simplest  example  of  a  fractal  is  the  so-called  "Koch  curve";  however 

the  derivation  of  D  =  log  4/log  3  =  1.2618  for  that  case  (as  well  as  the  discussion 

Afi 

of  many  other  examples  of  fractals)  will  not  be  given  here.  (See  Mandelbrot.  ) 

Thus,  fractals  are  of  fractional  dimension. 

Al  A2  A3 

It  has  been  shown  by  Lorenz,  Ruelle  et  al,  and  Feigenbaum,  that  if  a 

non-linear  system  has  an  attractor  in  its  phase  space,  which  is  also  a  fractal, 

then  the  system  can  exhibit  chaotic  motion.  It  follows  that  total  mixing  will  take 

place  in  coordinate  space  because  total  mixing  occurs  (in  the  fractal  sense)  in  the 

phase  space. 

The  following  citations  will  be  useful  for  further  reference.  In  particular, 

A 7  A8  A9  a  10 

the  works  by  Dold  et  al,  '  Swinney  et  al,  and  Treve  are  good  review 
papers  on  the  subject. 


A7.  Dold,  A.,  and  Eckmann,  B.  (1976/1977)  Turbulence  Seminar,  Springer 
Verlag,  N.  Y.  ,  156  pp. 

A8.  Dold,  A.  ,  and  Eckmann,  B.  (1976)  Turbulence  and  Navier  Stokes  Equation, 
Springer  Verlag,  N.  Y.  ,  194  pp. 

Af>.  Swinney,  H.  L.  ,  and  Gollub,  J.  P.  (1981)  Hydrodynamic  Instabilities  and  the 
Transition  to  Turbulence,  Springer  Verlag,  N.  Y.  ,  292  pp. 

All).  Treve,  Y.  M.  (l'>78)  Theory  of  chaotic  motion  with  application  to  controlled 
fusion  research.  Topics  in  Nonlinear  Dynamics,  S.  Jorna,  Ed,  , 
American  Institute  of  Physics,  pp.  147-220. 
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